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ABSTRACT 



Aims. We study the probability of close encounters between stars from a nuclear cluster and a massive black hole (10"* Mq < M. < 
10* Mq). The gravitational field of the system is dominated by the black hole in its sphere of influence. It is further modified by the 
cluster mean field (a spherical term) and a gaseous disc/torus (an axially symmetric term) causing a secular evolution of stellar orbits 
via Kozai oscillations. Intermittent phases of large eccentricity increase the chance that stars become damaged inside the tidal radius 
of the central hole. Such events can produce debris and lead to recurring episodes of enhanced accretion activity. 
Methods. We introduce an effective loss cone and associate it with tidal disruptions during the high-eccentricity phases of the Kozai 
cycle. By numerical integration of the trajectories forming the boundary of the loss cone we determine its shape and volume. We also 
include the effect of relativistic advance of pericentre. 

Results. The potential of the disc has the efffect of enlarging the loss cone and, therefore, the predicted number of tidally disrupted 
stars should grow by factor of 10^. On the other hand, the effect of the cluster mean potential together with the relativistic pericentre 
advance act against the eccentricity oscillations. In the end we expect the tidal disruption events to be approximately ten times more 
frequent in comparison with the model in which the three effects - the cluster mean field, the relativistic pericentre advance, and the 
Kozai mechanism - are all ignored. The competition of different influences suppresses the predicted star disruption rate as the black 
hole mass increases. Hence, the process under consideration is more important for intermediate-mass black holes, M. ^ \Q*Mq. 

Key words. Accretion, accretion discs - black hole physics - stellai' dynamics 



1. Introduction 

According to the present-day consensus, the observational ev- 
idence shows black holes populating the mass spectrum in two 
distinct intervals; stellar-mass black holes, 3>Mq < M, < 2QMq, 
that originate from collapse of massive stars and are revealed 
as members of binary systems; and supermassive black holes 
(SMBHs; 10*'Mo < M. < 10*^ Mq) hosted in nuclei of many 
galaxies including our own. The existence of intermediate-mass 
black holes (IMBHs; 10~Mq < M. < 10^ Mq) is an open is- 
sue, likewise their observational consequences and the form of 
the mass spectrum that might be filling the gap between the two 
well-established categories of black holes. 

Over the last four decades, a picture of galactic nuclei has 
emerged in which black holes are embedded in a dense stel- 
lar system and their masses are increasing by accretion of gas 
from an accretion disc, gradually inflowing towards the horizon 
(Begelman & Rees \WM Volonteri & Rees |2005l l. Remnants 
of stars tidally disrupted near a SMBH also contribute as a 
source of material (Zhao et al. 120021 ). Previously we demon- 
strated (Vokrouhlicky & Karas [T998l Subr et al. |2U04| ) that stars 
of the nuclear cluster can undergo episodes of very large orbital 
eccentricity if they interact with a self-gravitating disc near a 
SMBH. This mechanism sets stars on very elongated trajecto- 
ries with small periapses, although the overall tendency of the 
gas-assisted drag acts in the opposite way - towards the orbital 
circularization (Syer et al. 119911 Vokrouhlicky & Karas[T998]l. 



In this paper we study the fraction of stars that are set to 
highly eccentric orbits as a function of the black hole mass. We 
suggest this problem is relevant in the context of coexistence of 
a massive black hole with a surrounding cluster of stars. Exactly 
because of large elongation of the orbits, the fate of the rem- 
nant gas from the tidally disrupted stars is very uncertain and we 
do not attempt to solve this problem in its entirety. Instead, we 
concentrate our attention on a single aspect - the mutual gravi- 
tational interaction of stars and the gaseous disc. Our model is 
axially symmetric. We study how Kozai's phenomenon modi- 
fies and enlarges the black-hole loss cone and how this change 
depends on M,. We do take several subtleties of this scenario 
into consideration, namely the damping eff'ect of the star cluster 
and the relativistic pericentre advance of stellar orbits near the 
central black hole. These effects pose a potential threat to our 
mechanism. 

We find that a self-gravitating disc is capable of pushing 
more stars towards less massive (i.e. less-than-supermassive) 
black holes. Therefore, the effect should be particularly relevant 
for intermediate-mass black holes, provided they exist embed- 
ded within dense stellar systems and accrete from a gaseous disc 
or a torus (Madau & Rees l2005l Miller & Hamilton 2002) . 

There is a secondary motivation for our investigation; be- 
cause the efficiency of the mechanism discussed herein de- 
creases with the black hole mass increasing (it seems to be over- 
whelmed entirely by damping effects above M, ~ lO^M©), the 
present model may be particularly relevant for an ongoing de- 
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bate about the origin of ultra-luminous X-ray sources (ULXs) 
and their putative connection with IMBHs (for reviews, see van 
der Marel 2003; Fabbiano|2006ll. 

ULXs are extra-nuclear point-like X-ray sources with 
isotropic luminosities exceeding 10^^ ergs"'. It has been pro- 
posed (Colbert & Mus hotzky [19991 Makishi ma et al. l2U00l 
Matsumoto et al. 1200 U Portegies Zwart et al. 120041 Hopman 
et al. .20041 Baumgardt et al. 120061 ) that some of them may har- 
bour accreting IMBHs. Within the intermediate black-hole mass 
range the proposed mechanism is relatively efficient, and we 
can thus speculate that these ULXs might be triggered when a 
star is damaged near the tidal radius. This concept seems to be 
in line also with the idea that accreting pre-galactic seed holes 
may could indeed be detectable as ULXs (Madau & Rees 12005] ). 
Also, it appears to be in accord with the manifestation of a 
multi-colour disc component, as reported in spectra of a sam- 
ple of ULXs (Miller et al. l20D4t 120561 Fabian et al. l2004l l. and 
it also agrees with indications that ULXs are transient (Miniutti 
et al. 120061 ), possibly recurring phenomenon and that they are 
switched on during phases of active accretion (see Krolik l2004l ). 
However, the evidence for IMBHs is only circumstantial - nu- 
merous uncertainties persist in the mass estimates (see the dis- 
cussion in King et al. 120011 Gongalves & Soria 120061 ). and so 
our scheme is a mere speculation at this stage. 

2. Model and method 

2.1. Basic equations 

We assume that some kind of an accretion disc surrounds the 
central black hole and defines the equatorial plane of the system. 
The presence of the black hole provides a natural length-scale 
which we will use hereafter: Rg = GM,c^^ ^ 1.5 x 10^ M4 cm, 
M4 = M,/(10'*Mq). The disc is taken as axially symmet- 
ric and relatively light with respect to the central black hole 
(// = MilM, <*:!). Nevertheless, it is important that is greater 
than zero. To be specific, we consider two examples; (i) a con- 
stant surface density disc, characterized by its mass and the 
outer radius R^, and (ii) a limiting case of an infinitesimally nar- 
row ring of mass Mj and radius R^. The two cases serve as a 
useful test-bed because analytical expressions can be derived for 
the perturbation potential, Vi{R, z), and its first derivatives (Lass 
& Blitzer l 19831 Hure & Pierens |2U05T ). Notice also that the mas- 
sive ring surrounding the central black hole can represent a time- 
averaged system with a secondary black hole, sometime invoked 
in the framework of the hierarchical SMBH formation in cen- 
tres of merging protogalaxies (Volonteri et al. 120031 Giiltekin 
et al. 120041 ). Recently, Ivanov et al. ( 120051 ) applied the averag- 
ing technique to study tidal disruptions of stars in a galaxy cen- 
tre containing a supermassive binary black hole. However, these 
authors neglected relativistic effects, which exhibit growing im- 
portance near the black hole. 

The central cluster is bound gravitationally to the black hole 
and characterized by distribution function on the space of oscu- 
lating elements, Df{a, C\ , e, co), where C\ - rf cos ; is a normal- 
ized z-component of the angular momentum, 77 = (1 - e^)'^^; e 
is eccentricity, a is semimajor axis and o) is argument of pericen- 
tre. The orbit inclination / is measured with respect to the disc 
plane. We consider distribution of the form 



Df(a,Ci,e,w) = Ka^^'^ 



erj 



(1) 



with K being a normalization constant. Eq. ([T) represents the 
Bahcall & Wolf ( 119761 ) distribution with a random orientation in 
inclinations and linear distribution in eccentricities. The stellar 



cluster introduces a spherically symmetric perturbation to the 
central potential, Vc(r) oc r^l'^ . 

Length-scales of the model can be related to the mass ac- 
cording to the empirical relation M,(cr) (Tremaine et al. 120021 ) 
which we extrapolate in a naive way from galactic nuclei down 
to IMBHs range: M4 ^ \0'^<t\, where 0-2 = cr/(200km s"') is 
the velocity dispersion. This seems to be justified by kinematical 
studies of some globular clusters (e.g., M15, Gerssen et al. i 20021 
Gl, Gebhardt et al. l2005l 47 Tuc, McLaughlin et al. l2006l ) and 
it enables us to introduce the cusp radius, 

/?h = GM.o--^ = 2.25 X \{fM'^'^'\ = 0.1 1M{^^ pc . (2) 

We further assume the mass of the cluster M^iRh) - M,. Finally, 
we set the disc/ring radius to be equal to the characteristic radius 
of the cluster: R^ = Rh- Thenceforth only Md is a free parameter. 

The axial symmetry of the problem does not ensure conser- 
vation of the angular momentum; instead, only one of its com- 
ponents is conserved, so that the orbital eccentricity and inclina- 
tion of the cluster stars can secularly evolve and fluctuate. In our 
model we imagine that Kozai's phenomenorQ is responsible for 
oscillations of the orbital elements on the time-scale of 



Tk = 



A 



2nfi \ a 



1.6x 10-'M4. 



A /R, 



10^R„ 



yr. 



(3) 



The oscillation period exceeds the orbital period P(a) of each 
individual star, Tk » P = Ina^'^/iGM.y^^, by at least one 
order of magnitude. Factor A ^ I includes additional corrections 
arising from unaccounted effects. 

Precession due to an extended spherically symmetric poten- 
tial component accelerates the oscillations. For > O.IM, we 
find the correction term to the cycles of eccentricity oscillations. 

The period Tk from eqs. ([3])-(l4]) fits well with results of the nu- 
merical integration shown in Figure [1] where we evaluate the 
Kozai cycle period as a function of the mass of the spherical 
cluster. Tk in terms of orbital period P is determined numeri- 
cally by integrating the trajectories that reach the tidal radius R^ 
at the maximum of the eccentricity oscillations, i.e. Tk in Fig.[T] 
is given by the loss cone boundary trajectories. The dependency 
indicates that the shape of the loss cone is also modified, and we 
look to this more in the next section. 



2.2. Loss cone in fhe presence of Kozai's process 

The elongated trajectories bring stars of the nuclear cluster close 
to the centre. It is then natural to assume that passages near the 
black hole are moments important for the evolution of the system 
- feeding the black hole and triggering the accretion activity, 
as is commonly described in terms of loss-cone processes (e.g. 
Frank & Re es [T9761 Magorrian & Tremaine [T9991 Merritt & 
Wang l2005T ). 

A 'classical' loss-cone is determined by the process of tidal 
disruption of stars near the central black hole. Low angular mo- 
mentum trajectories are relevant because stars following these 
orbits have their periapses within the black hole tidal radius. 



1/3 



/?.-ioX"''( — 



-1/3 



Ri 



o 



R,, 



(5) 



' Originally (Kozai [79621 Lidov [T962t . the averaging method was 
applied to study the long-term evolution in the context of the restricted 
three-body problem. Naturally, the approach can be readily applied to 
discuss the motion in an axially symmetric potential. 
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Fig. 1. Kozai period Tk (solid lines) as a function of mass dis- 
tributed in the spherical potential of the nuclear cluster Time is 
scaled with the corresponding orbital period P. Principal param- 
eters of the system are the black hole mass, M4 = 10, and the 
mass of the ring, - O.IM,. Trajectories were characterised 
by C\ - 10"^ and integrated numerically. Two different values 
of semi-major axis are given by each curve. Also plotted is the 
approximation to Tk with the value of A from eq. Q, producing 
the slope -1/2 (dotted lines). 



where and are the stellar mass and radius. The loss-cone 
is emptied on a time-scale of ^ P provided the perturbing pro- 
cesses are not strong enough to deflect stars from their orbits. 

The above-described situation becomes more complicated 
when the timescales are modified in the presence of Kozai's 
mechanism. An 'effective' loss cone can be now associated with 
the new period Tk > P{a) operating in the system. The loss cone 
is formed by trajectories in phase space plunging below r - Ri 
by Kozai's oscillations. This view of the influence of the process 
takes into consideration the fact that orbital elements evolve in a 
systematical manner under Kozai's mechanism. In other words, 
eccentricities are pumped up to large values for extended peri- 
ods of time, which makes the effect distinct from the stochas- 
tic nature of gravitational scattering. Furthermore, the loss cone 
geometry is more complex because of the coupling that exists 
between the mean eccentricity and inclination of the orbits. 

We characterize the influence of highly eccentric orbits on 
the star capture rate by the fraction of stars that exhibit pericentre 
distances below a certain threshold radius, < Rmin- We denote 
this fraction TiiRmm, a, C\) and write it in the form 



1 rv'-w r^" 
'J^2=—\ de dw 0(eniax-emin) £>f(a,Ci,e,t<j), (6) 

L>2 Jo Jo 

where © is the Heaviside step function, e^i„ si- Rmin/o, and 
emax(fl, Ci,e, to) is the maximum eccentricity reached on a given 
trajectory. The normalization factor is equal to 



nyJl-C] pin 

02(0, C\)s j de I da) Df{a,Ci,e,oj) . 

Jo Jo 

Analogically we define integrated quantities 
1 f' 

TiiRn^; a) = — — dCi TiiRn^; a, C,) Diia, d) 
Di(a) Jo 



(7) 



(8) 



and 



5"(i?mi„) 



da Ti(Rmm,a)Di(a), 



where 



Di(fl) 



Jo 



dCi D2(fl,Ci). 



(9) 



(10) 



The function 'F{Rmm) can be interpreted as volume of the loss- 
cone normalised to the phase space volume occupied by the stel- 
lar cluster. (The correspondence between f and the volume of 
the loss cone would be exact in the case of uniform distribution 
function Df.) 

The integrals (l6Tl-(fT0t can be carried out analytically in the 
Keplerian case. In the non-Keplerian case, however, the function 
Smaxia, C\,e, (jj) is given in an implicit form which complicates 
the evaluation of the integrals. Our approach is thus based on di- 
rect numerical integration of the equations of motion of suitably 
chosen trajectories, as described below. 

2.3. Integration of the orbits 



According to the averaging technique (Arnold |1989t Brower & 
Clemence |19611 l the mean motion is taken over orbital period P. 
This allows us to study the long-term evolution of the osculat- 
ing parameters, (e, w, a, i) in an integrable (Keplerian) potential 
with a small perturbation. The approach relies on the existence 
of a third integral of motion (in addition to the energy and z- 
component of the angular momentum) which represents a mean 
of the perturbing part of the Hamiltonian. The orbital trajecto- 
ries form a congruence of closed curves in the space of mean 
eccentricity vs. argument of pericentre. 

We use the averaging technique to find parameters of the 
orbits that reach large eccentricities. Then, for those highly 
eccentric orbits we integrate the equations of motion directly. 
Accelerations we take in the form 



GM. 

a —r 



GM. 



4(r-v)v ■ 



4GM. 



v(yd + Vc), 



(11) 



assuming that the leading Newtonian term dominates. The 
bracketed term is the PNl correction to the gravity of the central 
mass (see sec. VII of Damour et al. 119911 ), while Va and Vc rep- 
resent axially and spherically symmetric perturbations, respec- 
tively. 

Figure |2] shows a set of tracks of different orbits with their 
corresponding values of a and Ci kept fixed. Pericentre argu- 
ment CO stands as the polar angle, whereas p{e) = (1 - 77)'^^ is the 
radial coordinate; maximum eccentricity is then on the perime- 
ter of these graphs, Cmax = (1 - Cj)'^^ = 0.995. Two classes of 
orbits can be distinguished - those which librate around to - njl 
(or to - 'injl) and circulating orbits for which the whole range 
of (0, 2n) is allowed. The most eccentric orbits attain their max- 
imum eccentricity e^ax at cos to - O 'm the outer circulating re- 
gion, as is evident from the plots. Hence, in order to evaluate the 
fraction TiiRmin',^, Ci) we integrate the orbit evolution starting 
from(e, w) = (emin, n'/2). The result of integration determines a 
boundary of the region occupied by trajectories with the maxi- 
mum eccentricity > gmin- Topology of the contour plot is crucial 
for the orbit behaviour 

The shaded area shows an intersection of the loss-cone with 
the (fl, Ci) plane. In the example of Fig. |2]we set = 0.989 



4 



V. Karas & L. Subr: Enhanced activity of massive black holes by stellar capture. 



1 




-1 I ' ' ' 1 I ' ' ' 1 I ' ' ' 1 

-1 -0.5 0.5 1 -1 -0.5 0.5 1 -1 -0.5 0.5 1 

P COS CO p COS CO p COS CO 



Fig. 2. Structure of the orbits is shown in (p, ai) polar plot. To find these orbits, we assume equations of motion in the form (fTTT i. 
For the gravitational field we take different terms into account, so that the effect of perturbations becomes apparent. Always present 
is the leading term of the Keplerian field (which dominates) and the axisymmetric field Wd (which represents the perturbation). In 
the left panel, no other terms are considered. In the middle panel, the PNl term of the central field is added. In the right panel, both 
PNl of the centre and Wc of the cluster are considered. High-eccentricity orbits pass through the shaded regions. Common values 
of the parameters; Me(i?h) = Md = O.OIM., R^^Ri.^ IQ^R^, Ci = 0.1 and a = QARd. 



and highlight the corresponding contour. This defines the bound- 
ary of the shaded area within which every orbit reaches e > e^nin 
at some moment of its evolution. In other words, this area con- 
tains highly eccentric orbits and it defines the size and the 
shape of the loss cone. Having the loss-cone boundary defined, 
we integrate Df{a,Ci,e,a)) over the shaded region and obtain 
T^iiRmin', a,Ci). Finally, we evaluate this function on a grid in 
the (a, Ci) space and integrate it numerically to find fiiRmin', a) 
and TiR^n). 

The three panels exhibit different possible topologies of the 
tracks of orbits that arise due to the perturbing forces. In order 
to construct these plots, the central Keplerian potential was su- 
perposed with the perturbation Wd (left); with PNl coiTection 
and Wd terms added (middle); and with PNl correction and 
V(yd + V'c) terms taken into account (right). 

In a conservative system, which is what we consider in 
Fig.|2] star tracks remain attached to the contours. However, ad- 
ditional dissipative processes might allow adiabatic evolution of 
the orbital parameters, so that stars can slide gradually across 
the level surfaces on very long time-scales. This may cause a de- 
cay of the orbits and help to bring stars to the centre. We do not 
consider such effect here, but see Subr & Karas (I2005I I where 
Kozai's resonance mechanism and the hydrodynamical (dissipa- 
tive) drag are both taken into account. 



3. Results 

As mentioned above, Kozai's resonance mechanism is active 
when the Newtonian central potential is perturbed by an axisym- 
metric term - i.e. the disc potential in our case, however, the 
effect is known to be fragile with respect to various other pertur- 
bations that may influence the motion of stars. Namely, it is sup- 
pressed when the central field is non-Newtonian, leading to the 
precession of the orbits. Therefore we want to clarify whether 
Kozai's mechanism is still relevant for the long-term dynamics 
of nuclear stars even if relativistic pericentre advance and the 
effect of the nuclear cluster are taken into account. 



3. 1 . Fractional probabilities 

First we neglect all perturbing terms in the gravitational field 
and assume purely Keplerian motion. The maximum eccentricity 
along each trajectory is simply en,ax(fl, Ci, e, w) = e - const. 
Assuming the distribution ([T]i, fractional probabilities Ti can be 
then written in the explicit form; 



TliRmm, a, Ci) = 0(?7min - Ci) , 

1 - Ci 

where T/min = (1 -e^in)'^^ 



T\(Rmm\a) - 77,, 



9P f?-^ 
^"min min 



(12) 



(13) 



/'("min) — ^lA (i^i 



5/4 _ 5/4 
«max flmin 



— 10 ^^min ^max ■ 



(15) 



To derive the approximation ( fTSl ) we assumed /?niin ^ flmin ^ 

^max- 

Under an axially symmetric perturbation of the central po- 
tential the total angular momentum of orbits is no longer con- 
served, nevertheless, its z-component remains an integral of mo- 
tion. This results in oscillations of eccentricity. An upper es- 
timate of the corresponding probabilities can be derived from 
the loss-cone condition, L- < Lmax, which can be expressed as 
emax(fl,Ci,e,c^) = {\-Cf)-^l^. Then, 



TliRmin, fl, Cl) = 0(/7min - Ci) , 
TliRmin, a) = 277niin - '?^n " ^ (2i?min 

and 



1/2 



(16) 
(17) 



,1/2 



_ 10(2^min)'^^(am^x-%-t) 
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Fig. 3. This figure exhibits a clear distinction between the cases with (on the right) and without (on the left) the relativistic pericentre 
advance taken into account. We draw the fraction TiiRmm, a, C\) of stars from the ensemble ([U that plunge below /?niin = 1007?g; 
a and C\ are constants of motion. The central black hole mass M, - 3.6 IQ^Mq, the disc mass = QAM,. In the bottom panels 
we assumed the central gravitational field perturbed by a narrow ring. In the upper panels, the case of a constant-density disc was 
considered. The left column is for the Newtonian potential of the central source, while on the right side we employed the post- 
Newtonian (PNl) model. Value of Ti is encoded by a colour palette, ranging from zero to unity. Therefore, this figure demonstrates 
how sensitive the value of Ti is with respect to the pericentre precession and the eccentricity oscillations - the two effects which 
act against each other Notice a sharp drop of below a certain critical value of a = Ot (in the panels on the right side). This break 
emerges once the pericentre advance is taken into account; it does not occur in the left panels because these were computed in purely 
Newtonian gravity. 



The last term on the right-hand side gives the values two to 
three orders of magnitude larger than those which follow from 
eq. dTsl l. The estimate ( fTsT i is good provided the Kozai mecha- 
nism is the only one determining the orbital changes. However, 
this resonance mechanism may be diminished by other pertur- 
bations, so the actual value of TiRmm) should be somewhere in 
between the two estimates. We will now discuss the importance 
of damping effects. 

The relativistic advance of pericentre has an overall tendency 
of decreasing the maximum eccentricity over the Kozai cycle. 
Below a certain threshold value of a the maximum eccentricity 
is not sufficient to bring stars inside the tidal radius; the loss- 
cone condition is then L < Ln,ax- We estimate this terminal value 
as 



a]^'iRlRlR-^^^^l-^ (19) 



(see AppendixlAt. For a > at the maximum eccentricity allowed 
by the relativistic corrections is large enough, so that the peri- 
centre gets below Rmm- Hence, the fraction of tidally disrupted 
stars is determined solely by the estimate ( fTlT l for a < a^, while 
above the fraction increases abruptly and saturates almost at 
the value (fTTI i. 



3.2. Stars plunging below a threshold radius 

In figureOwe plot the fractional probability TiiRmm, a, C\), as 
we obtained it numerically. Comparison between different pan- 
els clearly demonstrates the reason for the concern that the peri- 
centre advance of stellar orbits might completely erase the effect 
of resonance. 

One can see that Kozai's mechanism indeed sets stars on ec- 
centric trajectories reaching small radii and that the efficiency 
of the process increases with a decreasing down to a certain 
limiting value. We considered both the Newtonian and the post- 
Newtonian models of the central field, so the effect of the rela- 
tivistic pericentre advance can be distinguished. We find that ;F2 
is raised by a factor of ^ 100 due to the Kozai mechanism and 
it drops at small semi-axis; the exact value of a where that break 
occurs depends on the adopted form of the perturbing potential. 
Comparison of the panels reveals that the disc-like source of the 
gravity competes more successfully with the relativistic effect 
than a narrow ring of the same mass. 

The mean potential of the extended cluster of stars also 
causes secular precession of the orbits but this time it alters the 
results in a different way than the relativistic pericentre advance 
by attenuating the Kozai oscillations in the whole range of semi- 
major axes. Nevertheless, it still allows a substantial fraction of 
orbits to reach high eccentricities. 

The fraction Ti {Rt, a) is shown in figure|4] We find the damp- 
ing effects to be weakened if the mass of the axisymmetric com- 
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Fig. 4. Fraction fi (Ri; a) of stars plunging below the black hole tidal radius. Stellar orbits are perturbed by a massive ring of radius 
Rd and mass Md = i-iM, (indicated with curves). The dashed line is for the PNl approximation of the central mass gravity field; the 
effect of the outer stellar cluster is neglected. Solid lines represent two cases in which the cluster is taken into account (they differ by 
the disc mass parameter ju). Dotted are two power-law lines corresponding to the lower and the upper analytical estimates, eqs. ( [13]) 
and ([TtIi. respectively. Vertical dashed lines indicate the values of at (given by eq. iT% with u = 0.01 and 0.1, respectively); these are 
critical semiaxis values below which Kozai's phenomenon is suppressed. For a > at the eccentricity oscillations persist. Parameter 
values assumed in this plot are; in the left panel M, = lO'^M©; in the right panel M, - 10^ M0. 



ponent is higher, as expected. By increasing the mass parameter 
fi the curve of TiiRi, a) gets gradually shifted to larger values. 
Simultaneously, at is diminished. As consequence of these de- 
pendencies, more stars reach the tidal radius for lower mass of 
the central body than for higher masses. For the sake of definite- 
ness, we assume solar-type stars (M* - IMq and Ri, - IRq) in 
the evaluation of Ri in eq. (|5]). 

Finally, the overall fraction l^iRmin) is shown in figure |5] 
The limits on semi-major axes range are set mainly by numeri- 
cal arguments: the lower limit results from a sharp increase of the 
Kozai time when a is decreased. Below ami,, = O.OARh the or- 
bits would contribute to the tidal disruption rate only marginally. 
The upper limit has to be set due to limitations of our approach - 
at r ^ Ri the orbits exhibit chaotic motion and they no more fol- 
low the curves of figure|2] We assume that these orbits will lead 
to the disruptions and, therefore, results presented here should 
be considered as lower estimates. 

Line types help us to distinguish the importance of different 
effects in this plot. In particular: 

(i) Kozai's oscillations rise the value of 'F with respect to an 
unperturbed (fi - 0) case; 

(ii) PNl corrections to the central (Keplerian) potential reduce 
'F by factor ^2-^5 (the expected rate of tidal disruptions 
is diminished accordingly); 

(iii) precession due to self-gravity of the cluster decreases 'F 
further down by another factor ^10. 

As in the previous figures, it is clearly visible that the po- 
tential of the extended source, i.e. the disc, competes more suc- 
cessfully against damping than the ring. Increasing the mass of 
the axial perturbation not only rises F to higher values, but it 
also changes the slope. That tendency can be attributed to larger 
terminal value at of the semi-major axis, above which the rela- 
tivistic damping is negligible. 

The effect of relativistic pericentre advance increases with 
the mass of the black hole. The right panels of Fig.|5]demonstrate 
that F is anti-correlated with M, and grows towards less massive 
black holes (i.e. when going from SMBH to IMBH). 



4. Discussion 

We examined the idea of Kozai's effect assisted by an accre- 
tion disc as the mechanism enhancing the disruption rate of 
stars by a central black hole. This process is more efficient 
for intermediate-mass black holes than for supermassive ones. 
In our calculations we used PNl approximation of the central 
field. We performed these calculations also within the pseudo- 
Newtonian (Paczynski-Wiita 1 1 9801 ) framework (as a check of 
the orbits very near horizon) with very similar conclusions. 

The results for TiRt) can be directly interpreted as fraction 
of stars populating the loss-cone when an axisymmetric pertur- 
bation is applied on the cluster. We suggest that this picture is 
relevant in a system where gas clouds form an embedded disc- 
like structure with a non-negligible total mass near the sphere of 
dominance of a central black hole. 

4. 1 . Disturbing processes 

The view presented here is valid as far as other processes, act- 
ing on time scales shorter than Tk, do not manage to expell the 
stars from the effective loss cone. This would inhibit the Kozai 
process. One of the generic processes that may interfere with 
the influence of the axisymmetric perturbation is the two-body 
relaxation within the stellar cluster. 

Gravitational relaxation brings the diffusion time-scale into 
the problem, Tg = O^^T^, where 

oic = - ^Jrm « 1 (20) 

is an angular extent of the loss cone and 



' G^ClnAM^i* ' 

Here, «« is the number density of the stellar system. In A is 
Coulomb logarithm and C is a constant (in usual notation, 
ClnA ~ 10^; Spitzer 119871 ). In a cusp described by the 
Bahcall-Wolf distribution, velocity dispersion is comparable 
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Fig. 5. Fraction 'F(Rmm) of stars reaching the pericentre below a given radius r - /?niin- We assumed perturbing potentials Vd{R, z) of 
the disc with constant surface density (top panels) and a narrow ring (bottom panels). Left: the functional dependence of probability 
T'(Rmm) (soHd lines). Parameters compatible with Sgr A* were assumed: M. = 3.5 x IO^Mq, 0.01 < /U < 0.1, Rh - Rd - 
1.2 X 10^ Rg - 2 pc, and 0.04/?h < a < 0.4/?h. Right: the mass dependence, T(Rt(M,)), computed according to eqs. Q and (|5]l. In all 
panels dotted lines correspond to the analytical estimates ( fT4l ) and ( fTSl ); values given by the latter correspond within a factor of < 2 
to the case omitting both the relativistic pericentre advance and the stellar cluster reducing of the Kozai oscillations. Dashed lines 
represent models with only the relativistic precession switched on {/i - 0.01), while the solid lines include also the precession from 
the cluster potential (two cases are shown, /i = 0.01 and 0.1). 



with Keplerian velocity and the number density of stars de- 
creases with radius as «*(r) oc r"^'"*. Characterizing the number 
density by the stellar mass enclosed within the radius Rh, we 
get 



Tr 



10*'M: 



MX (MA I Rk 



5/4 



lO^R, 



1/4 



yr . (22) 



The effective loss cone which we have invoked in Sec. 12.21 
generally exceeds Gi^. From Fig.|5]we see that its volume TiRt) 
is (for // = 0.1 and a ring source of the perturbation) by factor 
> 10 larger than it is in an unperturbed system. Hence, appro- 
priate must be by a factor > 3 larger than Ox^. However, the 
geometry of the effective loss cone is more complicated. In par- 
ticular, it is narrower aiio - n/l where it actually coincides with 
the boundary of the classical loss cone. Therefore, at cj corre- 
sponding to the minimum eccentricity the new loss cone must 
be sufficiently wider so that its volume can exceed by a factor 
of ten the classical one. We estimate Sf^ ^ lO^k, which is also 
in accordance with our experience that the pericentre value on 



the boundary orbit of the effective loss cone typically changes 
by more than two orders of magnitude. The time scale on which 
stars diffuse across the effective loss-cone is then 



r' ^ lOOTfl ^ 0.04 M: 



\M.j \Moj \Rh 



-9/4 



(23) 



By comparing Tk and T'^ we deduce that for M4 > 10 the 
Kozai mechanism affects the angular momentum on a signifi- 
cantly shorter time-scale and we may consider the effective loss 
cone to be emptied by this process. 

For lower BH masses the interplay of Kozai mechanism 
and the two-body relaxation is more complicated, however, we 
still expect that systematic character of the angular momentum 
changes due to the Kozai mechanism will enhance the rate of 
tidal disruptions. A definitive study of these systems is probably 
not possible without employing high precision A^-body integra- 
tors. 

We imagine that consumption of stars by the central IMBHs 
is a transient process which enhances their activity during those 
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periods when stars are supplied into the loss-cone and then effi- 
ciently brought onto eccentric orbits. Permanent replenishment 
of this larger loss cone is not critical for the process considered 
here. Nevertheless, we may assume that the larger loss cone is 
replenished more frequently (also by other processes than the 
steady state relaxation - for example the effect of Brownian mo- 
tion of the central black hole, or a secondary black hole can con- 
tribute). However, details of the loss-cone refilling are beyond 
the scope of our present paper It will be interesting to com- 
pare these speculations with the results of A^-body simulations 
that have been recently applied to explore the tidal processes 
and their effect on IMBHs formation and feeding in dense star 
clusters (see Baumgardt et al. |2006t . 

To achieve a more complete description one has to invoke 
other diffusive processes acting together with the gravitational 
relaxation (e.g. the hydrodynamical drag by the disc gas that 
continuously helps re-filling the loss-cone; see e.g. Karas & Subr 
1200 ll for relevant time-scales) but, again, this is beyond the scope 
of our paper. 

4.2. Dependence on the black hole mass 

Let us consider two types of object which differ from each 
other by the black hole mass - (i) a supermassive black hole in 
the Galaxy center, and (ii) hypothetical intermediate-mass black 
holes that might reside in cores of dense star clusters. 

Precise tracking of the proper motion of individual stars is 
currently possible within an arcsecond area around the Galaxy 
centre (Genzel et al. 120031 Ghez et al. 12003b . and so the idea of 
applying our calculation to Sgr A* is quite natural and it raises 
a question about the origin of various perturbations that may act 
on stellar motion. Currently, ^ O.IM, can be set as an upper 
estimate of the mass of the molecular circumnuclear disc which 
extends from ~ 1 .5 pc to 3 H- 4 pc (Christopher et al. 12005 1 . 
i.e. on the outer edge of the black hole's sphere of influence. 
This is compatible with a narrow ring of radius = Rh - 2pc, 
and so we can directly evaluate the expected number of tidal 
disruptions. From Fig.|5]we read that the fraction ^ 5 x 10"^ of 
the total number of stars in the cluster with a,^^^ < 0.4/?d exhibit 
sufficiently large eccentricity oscillations. This corresponds to 
the mass in stars about ^ O.IM. ^ 3 x 1O''M0, thence we obtain 
^150 tidal disruption events per Tk - 10 Myr. 

So far we assumed that the source of the gravitational per- 
turbation was a gaseous disc or a torus, however, its origin could 
be different. In the Sgr A* there is a disc of young stars (Levin 
& Beloborodov l2003l Paumard et al. l2006l l and, possibly, a rem- 
nant gas from which these stars had been born. To model their 
effect we consider a thin disc of mass O.OIM, and surface den- 
sity oc R ~ extending between 0.03 ^ 0.3 pc. We found that ^ 2% 
of stars from the region < 0.3pc can undergo the oscillations 
that bring them down to R^. This means that during the period 
of Tk ^ IMyr after the formation of the gaseous/stellar disc 
there may have occurred up to ^ 100 more tidally disrupted stars 
compared to the result of calculations neglecting the effect. The 
phase of enhanced disruption rate may become prolonged if the 
orientation of the disc varies in time. This could be caused by 
precession in the outer galactic potential. 

In the case of IMBHs the resonance mechanism comple- 
ments the role of stellar encounters in the process of feeding the 
black hole. If the total mass of the cluster is about M, ^ lO'^M©, 
we expect the enhancement of the disruption rate ^ IMq per 
10^ yrs (with fj. = 0.1). Disrupted stars then provide material for 



accretion, thereby triggering a transient luminous phase of the 
object. 

We remark that the model has only formal validity near the 
lower end of the black-hole mass range because in this case the 
accretion rate of ^ lO^^'M, yr"' is comparable with the mass car- 
ried by individual stars; other processes must be important un- 
der such circumstances (which is consistent with the expectation 
that the accretion of stellar material onto IMBH is a non-steady 
process). Therefore, the scenario outlined above conforms to the 
assumption that ULXs are a transient phenomenon for which ac- 
cretion of the surrounding medium is essential. 

5. Conclusions 

Disruption of stellar bodies and subsequent accretion of the rem- 
nant gas are among likely mechanisms feeding black holes that 
are embedded in a dense cluster. We discussed one of the chan- 
nels that may contribute to this process. Kozai's mechanism can 
enhance the rate of such events, trigger the episodic gas supply 
onto the black hole, and, consequently, strengthen the activity of 
the system by raising the accretion rate. The process acts at char- 
acteristic time-scale of the Kozai cycle - typically Tk - 10^ yrs 
for IMBHs, during which the loss cone is depleted. 

Stars on highly elongated orbits are susceptible to tidal dis- 
ruption and hence they provide a natural source of material to 
replenish the inner disc. We have demonstrated that this phe- 
nomenon operates in the system even if it is disturbed by the 
relativistic pericentre precession and gravity of the nuclear clus- 
ter Therefore we can conclude that the presence of a gaseous 
disc of a small but non-zero mass, < < M,, helps dragging 
stars to the black hole, thereby feeding the centre and simulta- 
neously providing material that sustains and replenishes the disc 
itself. 
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Appendix A: The role of GR pericentre advance 

Here we want estimate the region of the parameter space where 
the general relativistic pericentre advance dominates over the 
Kozai effect. The quadrupole approximation leads to evolution- 
ary equations (Kiseleva et al. |1998l Blaes et al. 120021 ): 

Tk /? — = — ^ sin lai sin ; cos / , (A. 1) 

dr 8 

de 15 T y 
Tj^n — = — eT]~ sin2w sin / , (A. 2) 

df 8 

^-^'i-l W - 5 sin^ - [e' - sin^ 'J} - ^ | ■ (A.3) 

The last term on the right-hand side of eq. (IA.3b includes the 
relativistic pericentre advance. 

Equations ( lA.lb -l lATJt imply conservation of 

Ci-rj cos / (A.4) 

and 

Ci = (5 sin^ i sin^ w - 2) e" - . (A.5) 
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Solutions of two types can be found - circulating vs. librating - 
depending on the motion constants a,Ci, and €2- The circulation 
region exists always while the libration region occurs typically 
for small values of Ci . The two categories are separated by the 
separatrix curve that passes through e = point in the plane of 
(e, oS) polar coordinates. This can be used to determine the value 
of C2 on the separatrix: 



C2,s = - 



(A.6) 



Maximum eccentricity on the separatrix is given by equation 



which for rjs I simplifies to the form 



^r], — % 



(A.7) 



(A.8) 



We are interested in Ci small. For Ci < 0.01, the circulating 
solutions cover about half of the (e, w) parameter space. These 
orbits acquire eccentricities larger than at some moment dur- 
ing the orbit evolution. Hence, when e^ia, Ci) > 1 - a/R^in the 
fraction TiiRmm, o,C\) is close to the corresponding estimate 
( fT6] l which assumes that all orbits with given values of a and Ci 
reach the pericentre below 7? min. Setting 77^ - l-e^j^ - IR^in/a, 
eq. (lA.Sb gives an implicit formula for the boundary of these two 
regions in (a, Ci) plane. We can solve dA.St for Ci = to deter- 
mine the terminal value ( fT9] l of semi-major axis Of Below this 
value the relativistic pericentre advance inhibits the Kozai mech- 
anism. 
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